Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial.

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Genomic insights of evolutionary divergence and life history innovations in Antarctic brittle stars.

Understanding the drivers of evolutionary innovation provides a crucial perspective of how evolutionary processes unfold across taxa and ecological systems. It has been hypothesised that the Southern Ocean provided ecological opportunities for novelty in the past. However, the drivers of innovation are challenging to pinpoint as the evolutionary genetics of Southern Ocean fauna are influenced by Quaternary glacial-interglacial cycles, oceanic currents and species ecology. Here we examined the genome-wide single nucleotide polymorphisms of the Southern Ocean brittle stars Ophionotus victoriae (five arms, broadcaster) and O. hexactis (six arms, brooder). We found that O. victoriae and O. hexactis are closely-related species with interspecific gene flow. During the late Pleistocene, O. victoriae likely persisted in a connected deep water refugium and in situ refugia on the Antarctic continental shelf and around Antarctic islands; O. hexactis persisted exclusively within in situ island refugia. Within O. victoriae, contemporary gene flow linking to the Antarctic Circumpolar Current, regional gyres and other local oceanographic regimes was observed. Gene flow connecting West and East Antarctic islands near the Polar Front was also detected in O. hexactis. A strong association was detected between outlier loci and salinity in O. hexactis. Both O. victoriae and O. hexactis are associated with genome-wide increase in alleles at intermediate-frequencies; the alleles associated with this peak appear to be species specific, and these intermediate-frequency variants are far more excessive in O. hexactis. We hypothesise that the peak in alleles at intermediate frequencies could be related to adaptation in the recent past, linked to evolutionary innovations of increase in arm number and a switch to brooding from broadcasting, in O. hexactis.

Evolutionary innovations in Antarctic brittle stars linked to glacial refugia.

The drivers behind evolutionary innovations such as contrasting life histories and morphological change are central questions of evolutionary biology. However, the environmental and ecological contexts linked to evolutionary innovations are generally unclear. During the Pleistocene glacial cycles, grounded ice sheets expanded across the Southern Ocean continental shelf. Limited ice-free areas remained, and fauna were isolated from other refugial populations. Survival in Southern Ocean refugia could present opportunities for ecological adaptation and evolutionary innovation. Here, we reconstructed the phylogeographic patterns of circum-Antarctic brittle stars Ophionotus victoriae and O. hexactis with contrasting life histories (broadcasting vs brooding) and morphology (5 vs 6 arms). We examined the evolutionary relationship between the two species using cytochrome c oxidase subunit I (COI) data. COI data suggested that O. victoriae is a single species (rather than a species complex) and is closely related to O. hexactis (a separate species). Since their recent divergence in the mid-Pleistocene, O. victoriae and O. hexactis likely persisted differently throughout glacial maxima, in deep-sea and Antarctic island refugia, respectively. Genetic connectivity, within and between the Antarctic continental shelf and islands, was also observed and could be linked to the Antarctic Circumpolar Current and local oceanographic regimes. Signatures of a probable seascape corridor linking connectivity between the Scotia Sea and Prydz Bay are also highlighted. We suggest that survival in Antarctic island refugia was associated with increase in arm number and a switch from broadcast spawning to brooding in O. hexactis, and propose that it could be linked to environmental changes (such as salinity) associated with intensified interglacial-glacial cycles.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region.

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Combining genotypic and phenotypic variation in a geospatial framework to identify sources of mussels in northern New Zealand.

The New Zealand green-lipped mussel aquaculture industry is largely dependent on the supply of young mussels that wash up on Ninety Mile Beach (so-called Kaitaia spat), which are collected and trucked to aquaculture farms. The locations of source populations of Kaitaia spat are unknown and this lack of knowledge represents a major problem because spat supply may be irregular. We combined genotypic (microsatellite) and phenotypic (shell geochemistry) data in a geospatial framework to determine if this new approach can help identify source populations of mussels collected from two spat-collecting and four non-spat-collecting sites further south. Genetic analyses resolved differentiated clusters (mostly three clusters), but no obvious source populations. Shell geochemistry analyses resolved six differentiated clusters, as did the combined genotypic and phenotypic data. Analyses revealed high levels of spatial and temporal variability in the geochemistry signal. Whilst we have not been able to identify the source site(s) of Kaitaia spat our analyses indicate that geospatial testing using combined genotypic and phenotypic data is a powerful approach. Next steps should employ analyses of single nucleotide polymorphism markers with shell geochemistry and in conjunction with high resolution physical oceanographic modelling to resolve the longstanding question of the origin of Kaitaia spat.

Climatic change drives dynamic source-sink relationships in marine species with high dispersal potential.

While there is now strong evidence that many factors can shape dispersal, the mechanisms influencing connectivity patterns are species-specific and remain largely unknown for many species with a high dispersal potential. The rock lobsters Jasus tristani and Jasus paulensis have a long pelagic larval duration (up to 20 months) and inhabit seamounts and islands in the southern Atlantic and Indian Oceans, respectively. We used a multidisciplinary approach to assess the genetic relationships between J. tristani and J. paulensis, investigate historic and contemporary gene flow, and inform fisheries management. Using 17,256 neutral single nucleotide polymorphisms we found low but significant genetic differentiation. We show that patterns of connectivity changed over time in accordance with climatic fluctuations. Historic migration estimates showed stronger connectivity from the Indian to the Atlantic Ocean (influenced by the Agulhas Leakage). In contrast, the individual-based model coupled with contemporary migration estimates inferred from genetic data showed stronger inter-ocean connectivity in the opposite direction from the Atlantic to the Indian Ocean driven by the Subtropical Front. We suggest that the J. tristani and J. paulensis historical distribution might have extended further north (when water temperatures were lower) resulting in larval dispersal between the ocean basis being more influenced by the Agulhas Leakage than the Subtropical Front. As water temperatures in the region increase in accordance with anthropogenic climate change, a southern shift in the distribution range of J. tristani and J. paulensis could further reduce larval transport from the Indian to the Atlantic Ocean, adding complexity to fisheries management.

Global drivers of recent diversification in a marine species complex.

Investigating historical gene flow in species complexes can indicate how environmental and reproductive barriers shape genome divergence during speciation. The processes influencing species diversification under environmental change remain one of the central focal points of evolutionary biology, particularly for marine organisms with high dispersal potential. We investigated genome-wide divergence, introgression patterns and inferred demographic history between species pairs of all six extant rock lobster species (Jasus spp.), which have a long larval duration of up to two years and have populated continental shelf and seamount habitats around the globe at approximately 40o S. Genetic differentiation patterns reflected geographic isolation and the environment (i.e. habitat structure). Eastern Pacific species (J. caveorum and J. frontalis) were geographically more distant and genetically more differentiated from the remaining four species. Species associated with continental shelf habitats shared a common ancestry, but are geographically distant from one another. Similarly, species associated with island/seamount habitats in the Atlantic and Indian Oceans shared a common ancestry, but are also geographically distant. Benthic temperature was the environmental variable that explained most of the genetic differentiation (FST ), while controlling for the effects of geographic distance. Eastern Pacific species retained a signal of strict isolation following ancient migration, whereas species pairs from Australia and Africa, and seamounts in the Indian and Atlantic oceans, included events of introgression after secondary contact. Our results reveal important effects of habitat and demographic processes on the recent divergence of species within the genus Jasus, providing one of the first empirical studies of genome-wide drivers of diversification that incorporates all extant species in a marine genus with long pelagic larval duration.

Fine-scale population structure and evidence for local adaptation in Australian giant black tiger shrimp (Penaeus monodon) using SNP analysis.

BACKGROUND: Restrictions to gene flow, genetic drift, and divergent selection associated with different environments are significant drivers of genetic differentiation. The black tiger shrimp (Penaeus monodon), is widely distributed throughout the Indian and Pacific Oceans including along the western, northern and eastern coastline of Australia, where it is an important aquaculture and fishery species. Understanding the genetic structure and the influence of environmental factors leading to adaptive differences among populations of this species is important for farm genetic improvement programs and sustainable fisheries management. RESULTS: Based on 278 individuals obtained from seven geographically disparate Australian locations, 10,624 high-quality SNP loci were used to characterize genetic diversity, population structure, genetic connectivity, and adaptive divergence. Significant population structure and differentiation were revealed among wild populations (average FST = 0.001-0.107; p <  0.05). Eighty-nine putatively outlier SNPs were identified to be potentially associated with environmental variables by using both population differentiation (BayeScan and PCAdapt) and environmental association (redundancy analysis and latent factor mixed model) analysis methods. Clear population structure with similar spatial patterns were observed in both neutral and outlier markers with three genetically distinct groups identified (north Queensland, Northern Territory, and Western Australia). Redundancy, partial redundancy, and multiple regression on distance matrices analyses revealed that both geographical distance and environmental factors interact to generate the structure observed across Australian P. monodon populations. CONCLUSION: This study provides new insights on genetic population structure of Australian P. monodon in the face of environmental changes, which can be used to advance sustainable fisheries management and aquaculture breeding programs.

Mitochondrial and nuclear genetic analyses of the tropical black-lip rock oyster (Saccostrea echinata) reveals population subdivision and informs sustainable aquaculture development.

BACKGROUND: The black-lip rock oyster (Saccostrea echinata) has considerable potential for aquaculture throughout the tropics. Previous attempts to farm S. echinata failed due to an insufficient supply of wild spat; however, the prospect of hatchery-based aquaculture has stimulated renewed interest, and small-scale farming is underway across northern Australia and in New Caledonia. The absence of knowledge surrounding the population genetic structure of this species has raised concerns about the genetic impacts of this emerging aquaculture industry. This study is the first to examine population genetics of S. echinata and employs both mitochondrial cytochrome c oxidase subunit I gene (COI) and single nucleotide polymorphism (SNP) markers. RESULTS: The mitochondrial COI data set included 273 sequences of 594 base pair length, which comprised 74 haplotypes. The SNP data set included 27,887 filtered SNPs for 272 oysters and of these 31 SNPs were identified as candidate adaptive loci. Data from the mitochondrial COI analyses, supports a broad tropical Indo-Pacific distribution of S. echinata, and showed high haplotype and nucleotide diversities (0.887-1.000 and 0.005-0.008, respectively). Mitochondrial COI analyses also revealed a 'star-like' haplotype network, and significant and negative neutrality tests (Tajima's D = - 2.030, Fu's Fs = - 25.638, P < 0.001) support a recent population expansion after a bottleneck. The SNP analyses showed significant levels of population subdivision and four genetic clusters were identified: (1) the Noumea (New Caledonia) sample location; (2) the Bowen (north Queensland, Australia) sample location, and remaining sample locations in the Northern Territory, Australia (n = 8) were differentiated into two genetic clusters. These occurred at either side of the Wessel Islands and were termed (3) 'west' and (4) 'east' clusters, and two migrant individuals were detected between them. The SNP data showed a significant positive correlation between genetic and geographic distance (Mantel test, P < 0.001, R2 = 0.798) and supported isolation by distance. Three candidate adaptive SNPs were identified as occurring within known genes and gene ontology was well described for the sex peptide receptor gene. CONCLUSIONS: Data supports the existence of genetically distinct populations of S. echinata, suggesting that management of wild and farmed stocks should be based upon multiple management units. This research has made information on population genetic structure and connectivity available for a new aquaculture species.

Emerging patterns of genetic variation in the New Zealand endemic scallop Pecten novaezelandiae.

Both historical and contemporary processes influence the genetic structure of species, but the relative roles of such processes are still difficult to access. Population genetic studies of species with recent evolutionary histories such as the New Zealand endemic scallop Pecten novaezelandiae (<1 Ma) permit testing of the effects of recent processes affecting gene flow and shaping genetic structure. In addition, studies encompassing the entire distributional range of species can provide insight into colonization processes. Analyses of genetic variation in P. novaezelandiae (952 individuals from 14 locations, genotyped at 10 microsatellite loci) revealed a weak but significant regional structure across the distributional range of the species, as well as latitudinal gradients of genetic diversity and differentiation: estimates of migration rates supported these patterns. Our results suggest that the observed genetic structure and latitudinal gradients reflect a stepping-stone model of colonization (north to south) and emerging divergence of populations as a result of ongoing limitations to gene flow and insufficient time to reach migration-drift equilibrium. The low levels of interpopulation and interregional genetic differentiation detected over hundreds of kilometres reflect the recent evolutionary history of P. novaezelandiae and stand in contrast to patterns reported for other evolutionary older species at the same spatial scale. The outcomes of this study contribute to a better understanding of evolutionary processes influencing the genetic variation of species and provide vital information on the genetic structure of P. novaezelandiae.